Neutrino-induced fission of neutron-rich nuclei 
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We calculate neutrino-induced fission cross sections for selected nuclei with Z — 84 ~ 92. We 
show that these reactions populate the daughter nucleus at excitation energies where shell effects 
are significantly washed out, effectively reducing the fission barrier. If the r-process occurs in the 
presence of a strong neutrino fiuence, and electron neutrino average energies are sufficiently high, 
perhaps as a result of matter-enhanced neutrino flavor transformation, then neutrino-induced fission 
could lead to significant alteration in the r-process fiow in slow outflow scenarios. 



In this letter we calculate neutrino capture-induced fis- 
sion cross sections for heavy nuclei associated with the r- 
process. Matter-enhanced neutrino flavor transformation 
could enhance this process in the supernova/compact ob- 
ject environments commonly invoked as r-process sites. 
Recent observations of r-process abundances in low- 
metallicity, old galactical halo stars [1] show patterns 
which agree with the solar r-process distribution for nu- 
clides with mass numbers A > 130, but do not reproduce 
the solar r-process pattern for the lighter r-process ele- 
ments. In particular, these observed abundances show a 
peak around mass number A ~ 195, which follows the 
solar r-process distribution, and enhanced structures at 
around ^ 90 and ^ 132 which do not follow the solar 
pattern. 

It was recognized some time ago [2] that i/g capture 
on heavy nuclei in the post-collapse supernova envi- 
ronment would leave the daughter nuclei at the high 
excitation energies characteristic of Gamow- Teller reso- 
nances. This leaves these nuclei vulnerable to fission. Re- 
cently, Qian has demonstrated [3] that fission, induced by 
charged-current neutrino reactions within this neutrino- 
driven wind scenario [4], can account for the observed 
abundance patterns. In this model it is proposed that 
neutrino- induced fission occurs after the r-process freezes 
out (i.e. all initial neutrons are exhausted) and the pro- 
genitor nuclei decay to stability. It is further proposed 
[3] that no fission cycling occurs during the r-process, i.e. 
neutrino-induced reactions are unimportant during the 
r-process. 

Neutrino capture-induced fission cross sections have 
not been calculated before. Two aspects of nuclear 
physics conspire to make this process potentially impor- 
tant in dense environments with large neutrino fluxes: 

(1) the weak strength distribution in the charged current 
(neutrino capture) channel shows that the post-capture 
daughter nucleus will be left in a highly excited state; and 

(2) fission barriers are lower at higher excitation energy. 
It is expected that charged-current reactions on r- 

process nuclides will have larger partial fission cross sec- 
tions than neutral-current reactions, despite the fact that 
the latter can be induced by r neutrinos and their 



antiparticles, which, in a core bounce supernova explo- 
sion, might have larger average energies {{E^) ~ 20 — 25 
MeV) than the i^e neutrinos have {{E^} ^ 10 MeV). For 
the neutron-rich nuclei along the r-process path neutrino 
capture cross sections are quite large, as both allowed 
channels (Fermi and Gamow- Teller (GT)) are governed 
by sum rules which scale with the neutron excess N — Z 
and the Ve neutrino energy is large enough to excite the 
centroids of these allowed responses. Furthermore, the 
isobaric analogue state (IAS) and the GT centroid are 
located at energies in the daughter nucleus ~ 20 — 30 
MeV) which are significantly above the fission barriers in 
these nuclei. Hence, fission can represent an important, 
even dominant decay mode following neutrino-induced 
reactions on neutron rich nuclei. Obviously the principal 
competing decay mode is neutron emission, as neutron 
thresholds are also quite low in r-process nuclei. 

Our calculations of neutrino- induced reactions proceed 
through two steps. First we calculate the neutrino cross 
sections as functions of excitation energy in the final nu- 
cleus and then determine the decay mode of the final 
nuclear state using a statistical approach. The neutrino 
cross sections are calculated with the random phase ap- 
proximation (RPA), considering multipoles up to J = 4 
and both parities. (See refs. [5].) Our RPA scheme treats 
proton and neutron degrees of freedom separately and 
employs a partial occupancy formalism for non-closed 
shell nuclei. We adopt a zero-range Migdal force as a 
residual interaction. We note that the RPA satisfies the 
Fermi and Ikeda sum rules, which fix the total strength 
for the allowed transitions. 

In the second step we calculate for each final state with 
well-defined energy, angular momentum, and parity the 
branching ratios into the various decay channels using the 
stastical model code SMOKER [6], considering proton, 
neutron, a and 7 emission as well as fission. The fission 
barriers employed here were taken from the compilation 
of Howard and MoUer [7] and the neutron separation en- 
ergies from the mass table of Hilf et al. [8]. The final 
states in the residual nucleus were taken from the exper- 
imentally known levels supplemented at higher energies 
by an appropriate level density formula [6] . 
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Assuming a typical supernova v^, neutrino spectrum, 
i.e., a Ferm-Dirac spectrum with temperature = 4 
MeV and zero chemical potential, we have calculated 
the total (1^6,6") cross section and the neutrino-induced 
fission cross section for selected even-even nuclei with 
charge numbers Z = 84 — 92, covering the range from 
stability all the way to unstable r-process nuclides (Fig. 
1). As the total cross sections are dominated by allowed 
contributions (Fermi, GT), the cross sections increase 
linearly with neutron excess within the various isotope 
chains, simply reflecting the sum rules governing these 
two multipole transitions. Fission is an important decay 
mode, in particular for the Ra, Th, and U isotopes (Fig. 
1). The differences between the total and the fission cross 
sections, are mainly accounted for by the partial {v,e~n) 
cross sections, although for the lighter Po and Rn isotopes 
the decay into the gamma channel can compete with the 
fission decay branch. Due to the relatively high thresh- 
olds (Coulomb barriers), branchings into the proton and 
a channels are negligible. The competition between the 
two dominant decay modes, neutron emission and fission, 
are shown in Fig. 2 for selected Th and U isotopes. In 
our calculation, fission dominates the decay, except for 
the most neutron-rich nuclides shown. We note that this 
calculation only considers the decay branchings in the 
daughter nucleus, and does not follow multiple decays; 
i.e., it represents the "first-chance" fission cross sections 
[9]. 
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FIG. 1. Neutrino-induced charged-current cross sections 
on selected Po (upper panel), Rn (second panel), Ra (third 
panel), Th (fourth panel) and U (lower panel) nuclides. The 
total cross sections are shown by circles and the partial fission 
cross section by squares. A Fermi-Dirac spectrum with tem- 
perature = 4 MeV and zero chemical potential has been 
assumed for the Ve neutrinos. 



The competition between the dominant decay modes 
(fission, neutron emission) in a neutrino capture-excited 
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FIG. 2. Total (i/e,e~) (circles) and partial (i/e,e~n) (trian- 
gles) and neutrino-fission cross sections (squares) for selected 
Th (left panels) and U (right panels) nuclides. The calcula- 
tions have been performed for Fermi-Dirac neutrino spectra 
with temperature Tj, = 4 MeV and 8 MeV and zero chemical 
potential. 



daughter is governed by the relative values of the fission 
barrier Bf and the neutron separation energy Sn- The 
fission probability P/ is then approximately given by [17] 



1 + 4(m„/a^)i?2rexp {(% - 5„)/r} 



(1) 



where m„ is the nucleon mass and i? = 1.2 • A^l'^ is the 
nuclear radius. This formula assumes that the decaying 
nucleus is excited at energies E which are significantly 
larger than By and S'„. This is true in (i^ei e~) reactions 
with supernova neutrinos on heavy neutron-rich nuclei 
where i? ~ 25 MeV. Such excitation energies correspond 
to nuclear temperatures of T « 1 MeV in nuclei with 
A ^ 230 — 270. For simplicity we have assumed that 
T = 1 McV in the following. Eq. (1) yields P/ - 1/6 
if Bf ^ Sn, and P/ = 0.5 if Bf - Sn = -1.6 MeV. The 
difference U — B f — Sn is strongly dependent on the exci- 
tation energy. In fact, Eq. (1) is derived from statistical 
considerations involving the level density at vanishing nu- 
clear deformation (for the neutron emission probability) 
and at the saddle points of the double-humped fission 
barriers. The latter corresponds to a sizable nuclear de- 
formation, where the level density increases faster than 
at vanishing deformation. This reduces U with increasing 
excitation energy, enhancing the fission probability rela- 
tive to neutron emission [9,18]. This is consistent with 
the fact that the fission barrier in heavy nuclei is strongly 
influenced by shell effects [19,17], and these are washed 
out with increasing excitation energy. 

Using Eq. (1) we have inverted our calculated fission 
probabilities to obtain U{T) = Bf{T) — Sn{T), assum- 
ing T = 1 MeV. The desired quantity AU(T) = U{T = 
0) — U{T) is plotted in Fig. 3, where t/(0) has been de- 
rived from the compiled fission barriers [7] and neutron 
separation energies [8]. We note that the energy reduc- 
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tion is significant, amounting to about 4 McV on average. 
This result is in good agreement with earher estimates 
for heavy nuclei [9] . Although Fig. 3 shows some scatter 
among the studied nuclei caused by structure effects, we 
will assume that U{T) for supernova {i'e,e~) reactions 
on neutron-rich nuclei is lowered by 4 MeV compared to 
its ground state value. This allows for some interesting 
conclusions, which are rather independent of the chosen 
fission barriers and neutron separation energies. 
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FIG. 3. Tlie difference AU = U{0) - U{T) for neu- 
trino-induced fission by Ve neutrinos with a Fermi-Dirac dis- 
tribution of temperature = 4 MeV and zero chemical po- 
tential. U{T) = Bf{T) — Sn{T) is calculated from our partial 
cross sections using Eq. (1) and assuming a nuclear temper- 
ature T — 1 McV, while U{0) is derived from the tabulated 
fission barriers and neutron separation energies. 

For charged- current reactions with supernova neu- 
trinos one then has Pf = 0.5 for Bf - Sn ^ 2.4 MeV 
and Pf = 0.2 for Bf - Sn = 3.7 MeV, where Bf and S„ 
are the tabulated values appropriate for low excitation 
energies. Note that the predicted fission barrier heights 
vary quite significantly where modern evaluations (i.e. 
[20]) give higher barriers for neutron-rich nuclei than did 
earlier work (i.e. [7]). The recent fission barriers of [20] 
predict a fission probability Pf > 0.2 for most nuclei 
with Z > 92 and A > 230, with the exception of nuclei 
with lower Z-valucs around the potentially magic neu- 
tron number TV — 184. The fission barriers of Moller and 
Howard [7] allow for significant spontaneous fission prob- 
abilities for nuclei with Z > 87 and A > 230, including 
those around TV = 184. However, nuclei with A < 230 
have a smaller fission probability in neutrino- induced re- 
actions than tentatively assumed by Qian [3]. 

Nuclei on the r-process path have Sn ~ 1.5 — 2.5 
MeV. Such nuclei fission after excitation by neutrinos 
with probability Pf > 0.2 if Bf ~ 5.2 - 6.2 MeV. This 
condition is satisfied for some nuclei on the r-process 
path with Z > 96 for the fission barriers of [20] and 
for Z > 88 for those of [7], which are, however, likely 



too small for neutron- rich nuclei [6]. This would imply 
that neutrino-induced reactions can initiate a fission cy- 
cle, if the r-process production of superheavy elements 
occurs in a noticeable neutrino flucncc. Such a scenario 
might be conceivable, if neutrino oscillations occur. In 
contrast, the fission barriers of [20] are too high to allow 
for fission cycling during the r-process by /3-delayed or 
neutron-induced fission. 

We note that the typical fission cross section for Th 
and U isotopes (^ 400 • 10"''^ cm^) corresponds to a 
halfiife of ~ 0.08 s, assuming a neutrino reaction at a 
radius of 100 km above the neutron star and a typical 
supernova jy,, luminosity of 10^^ erg s"-'^. Such a halfiife 
is shorter than the expected halfiives for the r-process 
waiting point nuclei with N = 126,A> 195 [15] (and also 
with N = 184, A > 280 [16]). These typical halfiives are 
also comparable and may be shorter than the typical ~ 
0.1 s expansion timescale in "slow" neutrino-driven wind 
models. Thus, if the r-process occurs in a strong neutrino 
fiuence neutrino-induced fission on the progenitor nuclei 
during the decay to stability might affect the relative 
Th/U r-process abundance. This abimdance ratio is a 
necessary theoretical ingredient if one wants to deduce 
an age limit for the universe from the recently observed 
Th/U abundance ratios in old galactical halo stars [10]. 

The leverage that neutrino capture-induced fission has 
in an r-process set in a neutrino- driven wind is dependent 
on the i>e energy spectrum and on the neutrino fiuxes 
at the position where the neutrons are captured. Mod- 
els with an extremely fast outflow rate ( [11,12]) gener- 
ally have neutron capture occuring far from the neutron 
star where neutrino fluxes are low and, hence, neutrino 
capture-induced flssion effects could be scant, though 
post-processing fission could still be significant. 

Models with a slow outfiow rate suffer from a deficit 
of neutrons [2] associated with the "alpha effect." How- 
ever, these models can yield a viable r-process close to 
the neutron star if neutrino fiavor mixing effects are in- 
voked ( [13,14]). A hierarchical neutrino energy spec- 
trum, one where the mu and tau fiavor neutrinos are 
more energetic than the electron neutrinos remains a pos- 
sibility for at least some epochs following the bounce of 
the supernova core. In this case, matter-enhanced neu- 
trino fiavor transformation can play an important role 
in determining the efficacy of neutrino capture-induced 
fission, by making the average energies of the electron 
neutrinos larger and, hence, boosting fission probabili- 
ties in the region where neutrons are being captured in 
the r-process. (This is demonstrated in Fig. 2 for a T=8 
MeV neutrino spectrum.) Crudely, the relationship be- 
tween radial distance rg from the neutron star's center 
in units of 10 km and the temperature Tg in billions of 
Kelvins is « 22.5/(5iooT9), where Sioo is the entropy 
per baryon in units of hundreds of Boltzmann's constant. 
Typically, neutron capture in the "slow outflow" schemes 
takes place in the region where 1 < Tg < 3. The loca- 
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tion where a neutrino of energy Er will transform from 
muon/tau- flavor to electron flavor (and vice versa) is 

Tg'^sw « 1.3(20 McV/£;^)^/-\0.42/(ye + Y^))^/^ (2) 
X 5^0^0(3 X 10-3 eVV^m^ cos 26»)^^^ 

where (5m^ is the relevant difference of the squares of 
the vacuum neutrino mass eigenvalues (scaled here by 
the atmospheric neutrino value) and is the effective 
two-neutrino vacuum mixing angle, which for the z/g ^ 
i'^/t transformation channel in a strictly three-neutrino 
mass/mixing scheme would be roughly ^13 < 0.15. The 
experimental upper limit on this mixing angle precludes 
an adiabatic transformation at resonance in a straight 
MSW scheme, but we note that large effective matter- 
mixing could occur (depending on entropy and overall 
neutrino luminosity) on account of the flavor basis off- 
diagonal neutrino-neutrino forward scattering contribu- 
tions to the weak potential which determines neutrino ef- 
fective masses and matter mixing angles. In this equation 
Ye is the electron fraction and is the effective neutrino 
number fraction which enters into the neutrino forward 
scattering potential. Note that Y^, can be negative when 
the neutrino background is important. In fact, the neu- 
trino background potential can lead to near maximal fla- 
vor mixing in medium. Either way, the above expression 
and the implied location of significant fiavor transforma- 
tion is conservative: we may actually have a more ener- 
getic electron neutrino spectrum on account of " chaotic" 
maximal mixing well below this position. Crudely, the 
neutrino flux is (SxlG''^ cm-'^ s-'^r^r^ {lOMeV / {E^j)Lll. 
Here L^^ is the effective electron neutrino luminosity 
in units of 10^^ ergs s"^. If the entropy per baryon is 
5'ioo = 2, then the radius where a neutrino of energy 
Er = 20 MeV transforms is rg « 7 (corresponding to 
Tg w 1.6) and we would expect the typical lifetime 
against fission per big nucleus (inverse fission rate) to 
be Aj^^ ~ 0.05 s/L^^ where a = e,/i, r is the flavor of 
the progenitor of the electron neutrino when it leaves the 
neutrino sphere; whereas, if V,. +Y^ = 0.1, a possibility 
if neutrino mixing has been augmented by the neutrino 
background potential(s), then rg ~ 4.3 and Tg « 2.6 (for 
5ioo = 2) so that X^^ ~ 0.02 s/L^^. In either case, these 
lifetimes are shorter than typical waiting point r-process 
beta decay lifetimes and are shorter than at least a plau- 
sible range of expansion time scales, Tdyn ~ .015 s. This 
implies that the neutron capture flow could proceed out 
to some threshold nuclide mass in the 195 peak or just 
beyond, whereupon fission sets in, producing two fission 
fragments in the 130 peak, as outlined by Qian. 

To establish a steady state fission cycling scenario with 
neutrino capture-induced fission alone is problematic. If, 
in steady state flow, every seed nucleus is brought by neu- 
tron capture to a nuclear mass where the flssion cross sec- 
tion is greater than some threshold value, crj^, then flssion 
of this nucleus will result. Over a time At ~ 2rdyn there 



win be only some - 72 (Af/300 s^^) (At/0.03 s) (iV/8) 
neutrons liberated per threshold nuclear mass, where N 
is the assumed number of neutrons liberated per fission. 
Sustaining steady state fission cycle flow would require 
the liberation of some 70 to 100 neutrons per flssion frag- 
ment (mass ~ 130) and this is clearly untenable. Nev- 
ertheless, a more modest number of neutrons liberated 
per flssion coupled with the large rate of mass 130 flssion 
fragment production could represent a significant alter- 
ation in the r-process flow. At the very least it shows 
that the mass 130 and 195 peaks should have compara- 
ble abundances. 
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